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Abstract. We propose an original method for determining suitable refracting profiles 
between two media to solve two related problems: to produce a given wave front from 
a single point source after refraction at the refracting profile, and to focus a given 
wave front in a fixed point. These profiles are obtained as envelopes of specific families 
of Cartesian ovals. We study the singularities of these profiles and give a method to 
construct them from the data of the associated caustic. 
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1. Introduction 

Suppose that we are given a surface R of M 3 , which separates two media where the light 
travels with different velocities i>i and ti 2 , and a point F in the region corresponding to 
v\. In what follows, R will be called refracting profile and F point source of light. Let W 
be another surface in the region corresponding to v 2 , such that its normal lines coincide 
with the straight lines that emerge from F after refraction at R. The surfaces W with 
these properties are called wave fronts. A classical optical problem is to determine the 
wave front from the data of the point source F and the refracting profile R. It can be 
solved by using the Snell's law, or equivalently, the Fermat's principle. 

In this paper we consider the inverse problem, that is, given the wave front W and 
the source point F, we want to construct a refracting profile R such that W is the wave 
front obtained from a spherical wave front emanating from F after refraction at R. The 
same profile will solve the problem of focusing the wave front W in the point F. 

We will give a method based on the optical properties of the Cartesian ovals (section 
2), to obtain a family {-R a } of refracting profiles parametrized by a non-negative real 
number a, with the property of producing normal lines to W after refraction at each 
profile of the rays that emerge from F. Each profile R a will be the envelope of a family 
of Cartesian ovals (see figure 1). A similar method have been used in |XJ to construct 
reflecting profiles as envelopes of certain families of conies. 

Next section contains the details of the construction of the profiles R a and the proof 
that, under very general conditions, these profiles are immersed surfaces of M. 3 that have 
singularities only if they contain centres of curvature of the wave front W. Moreover, 
we prove that the singularities of the families of the profiles R a sweep out the caustic of 
the wave front W. In section 3 we study the optical physical sense of these profiles. 

Finally, in section 4 we show an interesting relation between the profiles R a and the 
caustic of the given wave front, which also allows us to construct them from the caustic. 

2. Construction of the refracting profiles 

Let us remember some properties of the Cartesian ovals. The ovals of Descartes or 
Cartesian ovals were introduced by Descartes in 1637 in his Dioptrique, dedicated to 
the study of light refraction. A Cartesian oval is the locus of the points from which the 
distances r\ and r 2 to two fixed points F\ and F 2 , called foci, verify the bipolar equation 

air 1 + a 2 r 2 = fc, (2.1) 

where a±, a 2 and k are constants. Observe that this equation includes the bifocal 
definition of conies for the particular cases of k > and a\ = a 2 > (ellipse), or 
a\ = —a 2 > (hyperbola). 

A Cartesian oval has a third focus F 3 , and the oval can be defined by any two of 
the foci (see j2]). In particular, when | ai/a 2 \ goes to 1, the third focus goes to infinity 
and we get the conies. 
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Figure 1. The wave front W focuses in F after refraction at the profile R a . This 
profile is the envelope of a family of interior Descartes ovals 0% with foci F and x £ W; 
ni and ri2 are the refractive indices. 



The so-called complete Cartesian oval is the set of curves associated to the bipolar 
equation 



Only two of the four equations obtained from these double signs are not empty. If 
a\ 7^ a 2 , they are closed curves, one interior to the other (see figure 2 with F\ — F and 
F 2 = x). These two curves intersect only when two foci coincide. In this case we get 
the so-called Limagon of Pascal. See |2j for a detailed study of these ovals. Reference 
[3] contains some applications of these ovals to holography. 

A complete Cartesian oval is symmetric with respect to the straight line determined 
by the foci F% and F 2 . Then, if we consider F\ and F 2 belonging to K 3 , the bipolar 
equation (|2.2|) is the equation of a surface of revolution in M. 3 with axis of revolution FiF 2 . 
Since the intersection of this surface with any plane containing F\ and F 2 is a complete 
Cartesian oval, we have called this surface complete Cartesian oval of revolution. 

The fact that makes interesting Cartesian ovals of revolution in the context of 
geometrical optics is the following. Consider a surface R separating two media with 
refractive indices n\ and n 2 , and two points F\ and F 2 in medium n\ and n 2 respectively. 
If we want that a ray coming from F\ goes to F 2 , then the point P where the ray cross 
S must verify that n\\F\P\ + n 2 \F 2 P\ be constant, but this shows that S is part of 
a Cartesian oval of revolution, as we see looking at Eqn. (j2.2j) . Because the refractive 
indices n\ and n 2 are related with the light velocities V\ and v 2 by rij = c/vi, i — 1,2, 
where c is the light vacuum velocity, the above condition is equivalent to say, that all 
the rays spend the same time going from F\ to F 2 , and this is the Fermat principle. 
We will see in Section 3, that only some relative position of the foci with respect to the 
complete oval of revolution and some portions of this oval have physical sense. 



ai7T ± a 2 r 2 = ±k (k > 0). 
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In what follows let O® denote the complete Cartesian oval of revolution with foci 
F and x £ W, satisfying Eqn. (j2.2j) for the parameters a\ and 02 equal to the refractive 
indices n\ and ri2 respectively, and the optical path length between F and x, k, equals 
to 2a. Let R a be the envelope of the family {0%.}xew- We will give an explicit 
parametrization of R a and from it we will determine which portion have physical sense 
as refracting profile. 

From now on we assume that F is the origin of R 3 , so any point y £ must verify: 
ni\y\ ± 77.2 1 2/ — x\ = ±2a. (2-3) 

Therefore, for any a > 0, we can define the family of complete Cartesian ovals 
of revolution {0®} x ew by mean of the maps f a : W x R 3 — > E given by f a (x,y) = 
n \\/V ' V =t n 2\/ (y — x) ■ (y — x) ± 2a, so that for any fixed x £ W the zero set of the 
map f a>x : R 3 -> R given by / ajl .(?/) = f a (x,y), is the Cartesian oval O". 

According to the definition of the envelope of a two-parameter family of surfaces 
(see Refs. |U E]), the envelope R a of the family {OJ} l6 ^ is the set of points y £ R 3 
such that exists t = (ti, t 2 ) £ £7 verifying the equations: 

fa{x{t x M): y) = 0, 7^(x(ti, t 2 ), y) = 0, -^(x(h, t 2 ), y) = (2.4) 

Oti ot 2 
for t = (ti,t 2 ) being the standard coordinates in R 2 and x : U C R 2 — > R 3 a 
parametrization of W. 

The envelope R a is therefore obtained by eliminating t\ and t 2 between the following 
equations 

fMt),y) =0 (2.5) 

Then, equation f|2.6|) becomes 

(y-x(t))- — = 0, i = 1,2, (2.7) 

which, is also the equation of the line r n (x), normal to W at x(t), and thus it follows 
that 

R a = (J °" nr n(x). (2.8) 

This proves that the normal line to W at x intersects the envelope R a at the points y(x) 
where the envelope itself is tangent to the Cartesian oval 0£ (see figure 2). 

In order to give an explicit parametrization of R a , we suppose that W is oriented 
and we take n(x) to be the outward unitary normal vector field to W. According to 
(|2.8)l . R a can be given by the parametric equation 

y = x + X(x)n(x) (2.9) 

where A(x) has to be determined by the condition that y £ 0°, that is, y has to verify 
(EL 
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Figure 2. The wave front W and the four profiles Rfj = 1,2), which are the 
envelopes of the interior (0° J and exterior (0% e ) Descartes ovals. These ovals have 
foci F and x € W; x' is the third focus. The singularities of the profile i?f 2 are on the 
caustic C of the wave front . 



A straightforward calculus gives four A's: 

x _ 2an 2 - (-l)^(x ■ n) + (-l)^y/A~ 

Ay— -2 — — 2 , i — -M, l^-iuj 

' 4 2 '•'1 

which are defined only when Aj = (2cm 2 — (— l) l n\(x ■ n)) 2 — {n\ — n 2 )(4a 2 — n\x 2 ) is 
positive. Then R a can be decomposed in four sheets R a = Ujj = i i2 -R"j, where R?j is 
parametrized by 

y = x + Xi ) j(x)n(x) (2-H) 

Figure 2 illustrates these four sheets in a 2-dimensional example. 

Let us now see that the map yij : U C M 2 — > M 3 given by 2/y(t) = + 
Ajj(x(f))n(x(f)), = 1,2) is an immersion, so that _R a = Uij^^-Ry is an immersed 
surface of M 3 . 



To simplify the notation we avoid the indices i,j in (|2.11jl . To see that y is an 
immersion we have to verify that the vectors 
dy = dx +x dn + d\^ k = 
dt k dt k dt k dt k 



are linearly independent. 

3 noromafri^ofiAn l m 1/1/ txro r»on QOCiTm£3 f hof f no ~t ronf atc .. 



If we chose a suitable parametrization in W, we can assume that the vectors ° J '' 



are orthonormals. With respect to the orthonormal base { J-f-, J|f-, of M 3 , the two first 
components of the vectors J^-, form a (2 x 2)- matrix given by: 

Sl ^ ^ +A |£.^ = ^_ A ^. n y (2 , 3) 



where for the second identity we have used that = f • nj — • J^- + gf^. • ft- 
This matrix is singular if and only if 1/A is an eigenvalue of the matrix ( ■ n^j which 
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is the matrix associated to the second fundamental form of W at x(t) with respect to the 
base {J-p-, J^-} of the tangent plane to W (see p. 141 of 6J). Therefore, its eigenvalues 
are the principal curvatures of W at x(t). 

It follows that the matrix ()2.13|) is singular if and only if 1/ A is a principal curvature 
of W at x(t), or equivalently, y(t) is a centre of principal curvature of W at x(t). Then 
we have proved the following theorem: 

Theorem 2.1. Suppose that R a does not contain any centre of principal curvature of 
W . Then R a is an immersed surface ofM 3 . 

According to the above theorem, the singularities of the refracting profiles R a are 
the points y G R a which are centres of principal curvature of W. The locus C, of the 
centres of principal curvature of the wave front W, or equivalently, the envelope of the 
normal lines to W, is called the caustic of W. 

In what follows we only consider wave fronts whose caustics are non-degenerate in 
the sense that each point c of the caustic corresponds to an unique point x of the wave 
front . 

It is known that when the wave front propagates its singularities slide along the 
caustic, see [HJ- It was proved in [I] that a similar result holds regarding the 
singularities of the reflecting profiles. In the next theorem we prove that an analogous 
result is valid for the profiles R a . 

Theorem 2.2. For each c EC ofW there are two refracting profiles R ai and R a2 such 
that c is a singularity of both of them. 

We express this fact saying that the singularities of the family {-R a } a >o sweep out 
the caustic C of W twice. 

Proof For each point c G C take a x = |rai|c| + ^ C ~ X W and a 2 = |ni|c| -"^H . □ 

Figure 3 gives an illustrative example of how the singularities of the profiles R a 
sweep out the caustic C. 

3. Physical sense of the different profiles 

We will study now which sheets of R a are adequate as refracting and focusing profiles. 
That is, which sheets produce the wave front W after refracting at them a spherical 
wave front emanating from F, or equivalently, focus W in F. First of all, we will give 
an explicit parametrization of the family of complete Cartesian ovals {O^xew- 

To reduce the problem to the 2-dimensional case, we consider for every x G W the 
plane n(F, x) determined by F, x, and n(x). Now if we take F as pole and Fx as polar 
axis, the polar equation r = r(x, tp) for a point y G 0% is determined by the condition: 

\n\T ± 712^1 = 2a (3-1) 

where r = \y\, r 2 = \y — x\, and n\ ^ n 2 . 




On the other hand, r and r 2 are, together with \x\ the sides of a triangle. Then: 

r\ = r 2 + \x\ 2 — 2r\x\ cosp. (3-2) 
Using (|3.1|) and (|3.!2j) we get two solutions: 



2an 1 -n 2 2 \x\cos(f± y/A(x,(p) 
r±{x, if) = - 2 — - 2 (3.3) 

where 

A(x, if) = (2ani — n\\x\ cosy?) 2 — (n\ — n\){Aa 2 — n\\x\ 2 ). (3.4) 

Consider now a cartesian coordinate system of the plane ir(F,x), F,X,Y, with 
origin at F. With respect to this coordinate system the cartesian equations associated 
to r±(x, ip) are: 

V± = {r±{x, if) cos(y9 + 6), r±(x, ip) sin(v? + 6)), (3.5) 

where 9 is the angle determined by x and FX. 

In figure 4(a) we have represented an example with these two solutions. It illustrates 
also the following discussion. 

The solution r + (dashed line in figure 4(a)), has two parts: one in the interior oval 
0£j, and the other in the exterior oval 0% e - They correspond to the values of <p for 
which A(x, <f) > 0, and are separated by other two regions where A(x, if) < 0. So there 
are four values of <f for which A(x, if) = 0, as can be deduced from ()3.4|) . 

These four values of if are the corresponding to the values of the angle between 
FX and the four tangent lines to the interior O a x i and the exterior O a x e ovals passing 
through F. 

For the solution r_ (solid line in figure 4(a)), we have an analogous decomposition. 
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Then both the interior and exterior ovals are made up of two pieces, one 
corresponding to r + (dashed line in figure 4(a)) and the other corresponding to r_ 
(solid line in figure 4 (a)). 

When F is in the interior of the ovals, A never vanishes, so in this case, the exterior 
and interior ovals coincide with the ovals r + (dashed line in figure 4(b)) and r_ (solid 
line in figure 4(b)) respectively. 

A sheet of R a could have physical sense as refracting or focusing profile only when 
F and W are at different sides of the tangent plane to the sheet at any point of it, and 
a point y being in that sheet must verify n\\y\ + n 2 \y — x\ = 2a. On the other hand, for 
any x G W, a point y in the interior oval 0%i verifies (see 

ni\y\ + ri2\y — x\ = 2a (3.6) 

and a point y of the exterior oval e verifies 

— rii\y\ + n 2 \y — x\ = 2a. (3.7) 

From condition ()3.6|) we know that the refracted ray at the point ?/ G 0"j of a ray 
emanating from F, passes through x (see figure 2). We know also by (2.7) that the 
normal line to W at x intersects the envelope _R" 2 at the points y(x) where the envelope 
itself is tangent to the oval 0%. This fact together with the above observation means 
that the rays produced from the point source F by refraction at R® 2 are the normal 
rays to W. Therefore we have proved the following result, which can be considered the 
main result of this paper. 

Theorem 3.1. Let W , F , and Rf 2 be as defined above. Then the profile Rf 2 refracts 
the rays emerging from F into the normal rays to W. 

In other words the theorem says that the wave front W can be obtained from a 
spherical wave front emerging from F after refraction at R\ 2) or equivalently the profile 
R12 focuses the wave front W in F. 

We can prove with an analogous reasoning that the sheet R 2 2 has a physical sense 
when we consider F as a virtual source point. Let us assume that the profile R 22 
separates two media with refractive indices n 2 in the side containing F and n\ in the 
other side, and that both F and W are at the same side of the tangent plane to the 
sheet for any point of it (see figure 4(a)). Then we have the following: 

Theorem 3.2. The wave front W , when propagates towards R 22 and after refraction 
at it, produces a spherical divergent wave front S with centre F. 

Remark. If we reverse the direction of propagation of the wave front, the above 
theorem can equivalently be stated as: A convergent to F spherical wave front S, when 
propagates toward R 22 , and after refraction at it, produces the wave front W. We say 
in this case that F is a virtual source of W for the refracting profile R 22 . 

Proof. Consider a point z G S, so it verifies \z\ = R (constant). From the refraction law 
the optical path length n 2 \y — x\ + n\\z — y\ from x G W to z G S must be constant. 
Moreover, since \z — y\ = R — \y\ then — n\\y\ + n 2 \y — x\ must be constant, and this 
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Figure 4. (a) If F is in the exterior of the ovals, a spherical wave front emerging from 
F produces, after refraction at 2> the wave front W, and, after a posterior refraction 
at i?2 2 , appear again the initial spherical wave front, (b) If F is in the interior of the 
ovals, a spherical wave front converging to F produces, after refraction at R® 1 , the 
wave front W, and a posterior refraction at i?f 2 produces a spherical divergent wave 
front with centre F. 



is just the equation that satisfy the points y in the exterior ovals {0® e } xe w, whose 
envelope is R 22 . ^ 

Using together theorems Kill and IH.2I we can derive an interesting optical property. 
Consider a wave front W, a source point F, and a parameter a such that W and F are 
at different sides of R® 2 . Suppose given refractive indices n\ < n 2 such that the region 
delimited by R\ 2 and R 2 2 is filled with a medium of refractive index n<i , and the exterior 
of this region with refractive index ri\. Then a spherical divergent wave front emerging 
from F produces, after refraction at Rf 2 the wave front W, and after another refraction 
at R 22 , produces a propagation of the initial spherical wave front (see figure 4(a)). This 
process is a special kind of the known as do-nothing machines (jHj), and it can have 
useful applications to optic, acoustic and radar devises. 

The physical sense of the different sheets of the profiles R a depends on the kind of 
ovals (exterior or interior), whose envelope is the sheet, as well as the relative position 
of the considered sheet with respect to the wave front W and the source point F . As 
has been shown before, when W and F are at different sides with respect to the profile 
Ri2, then R1 2 and R 22 have got the suitable properties of refracting profiles. 

There is another profile sheet with physical sense, namely the profile Rx V Consider 
a value of the parameter a such that F and W are in the same region with respect to 
the tangent plane to this sheet, i.e. n(x) ■ y > 0, for n(x) pointing to the medium rii 
(see figure 4(b)). Reasoning as before, we obtain that under these condition: 
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Theorem 3.3. Any spherical wave front S, convergent to F from the medium of n\, 
after refraction at R^ ± produces the wave front W . 

Equivalently, the theorem says that the wave front W after refraction at R\i, 
produces a spherical divergent wave front S with centre F. 

A nice optical property can be deduced from theorems 13.21 and 13.31 

Consider that the region delimited by Rf l and R 2 2 * s filled with a medium of 
refractive index n 2 , and the exterior of this region with refractive index n\. Then a 
spherical wave front converging to F and after refraction at Rfi, produces the wave 
front W, and subsequent refraction at R 22 gives a divergent spherical wave front with 
centre F (see figure 4(b)). This property can also have useful applications. 

Notice that if we want to have physical sense, in all the above discussion we must 
exclude the region of the profile in which there are total reflection, that is, we have to 
exclude the points y of the profile such that the corresponding incidence angle 9\ does not 
give a real value for the angle of refraction 9 2 . It occurs when light is propagated from 
an optically denser medium into one which is optically less dense, i.e. when 712 < n%, 
provided that the incidence angle 9\ exceeds the critical value 9 C given by sin 9 C = n%jnx. 

4. Construction of the _R a -profiles from the caustic 

In this section we develop a method to construct the i? a -profiles of a given wave front 
W with respect to a point F, starting from the data of its caustic C. As we set in the 
above section, we only consider wave fronts whose caustics are non degenerated in the 
sense that each point c of a caustic corresponds to a unique point x of the wave front. 

For each regular point c G C corresponding to a point x G W and a principal 
curvature radius p, such that c can be expressed as c = x + pn(x), set a' = a + |n 2 and 
a" — \a — I n 2 \. 

Consider the families of ovals of revolution {0"'} ce c and {0""} cg c, with foci F and 
c and parameters a' and a" respectively. 
Then we have: 

Theorem 4.1. With the notation above, the R a -profiles are in the envelopes of the 
families of complete ovals of revolution {O"' U ^c"}cec- 

Proof. We shall prove the result for the convex regions of the wave front. The proof for 
the concave regions is analogous. 

In the convex regions the sheets Ri X , R\2i ^21; anc ^ -^22 are ^ e envelopes of 
the families of ovals {0"' e } c6 c, {Oc'i}cec, {Oc" e }cec, an d {O^'jcec respectively; here 
the letters e and % denote exterior and interior sheet respectively of the corresponding 
complete oval. We will centre our attention to the case of refracting profiles, that is, 
Ri2- The proof is similar for the other cases. Figure 5 illustrates with an example the 
families of ovals whose envelopes are and Rf 2 - We have to prove that Rf 2 is the 
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Figure 5. (a) For the case of W convex, Rf 1 and R1 2 are the envelopes of the ovals 
0° e and respectively, which have foci F and c E C. (b) For a concave W, R1 1 
and Rf 2 are the envelopes of the ovals 4 and e respectively, which have foci F 
and c S C. 



envelope of {0"'j} c6 c. We already know that the refracting profile is the envelope of the 
interior ovals {0"j} cG iy, verifying 

rii\y\ + n 2 \y — x\ = 2a. (4.1) 

On the other hand the points y G verify 

ni\y\ + n 2 \y — c| = 2a' = 2a + pn 2 = 2a + \x — c\n 2 , (4.2) 

where c is the centre of curvature of W at x. We shall show that and are 
tangent at the point y in which is tangent to the refracting profile (see figure 5(a)). 

In order to demonstrate this result it suffices to show that the intersection between 
0"'j and 0°j is only one point, precisely the point y where 0° Ei intersects the normal 
line to W at x. This follows from the fact that y verifying both (|4.1|) and ()4.2|) implies 
\y — a; | = \y — c| — \x — c\, what it is possible only if y is in the straight line determined 
x and c. Observe that \y — c\ > \x — c\ since for a convex wave front the caustic C and 
F are in opposite sides with respect to W (see figure 5(a)). 

For the concave case, R% l7 R\2i R211 an d -^22 are the envelopes of the families of 
ovals {Oc'e}cee, {0"'-} ce c, {0<} ceC , and {0"' e } cgC respectively. Figure 5(b) illustrates 
this case for Rf 1 and R1 2 . □ 
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5. Conclusions 

Using envelopes of families of Cartesian ovals of revolution, we have given an original 
procedure to solve the following inverse problem: find a refracting profile R separating 
two media such that for a given wave front W and a fixed point source F, the wave front 
produced by F after refracting at R is precisely W . In fact, with the given procedure we 
have constructed a one-parameter family of profiles, each member of this family having 
four sheets. We have studied the optical physical sense of these sheets. 

We have proved that these refracting profiles are immersed surfaces in R 3 , with 
singularities in the points corresponding to the centres of curvature of the wave front. 
Finally, we have given another way to reconstruct the refracting profile from the caustic. 

The proposed method can be of interest in the study and design of refracting 
profiles. The fact that each one of these profiles is constructed as the envelope of 
Cartesian ovals of revolution could be used to design it as a tesselation of smaller ovals 
profiles. This composed profile could be a good approximation of the original one, and 
therefore could be useful in the construction of optical, radar, or acoustical devices. We 
have also found that with a suitable choice of two of these profiles, it is possible to 
construct devises with do-nothing machines properties. 

For ri\ = n 2 the Cartesian ovals are ellipses or hyperbolas, and some of the sheets of 
the i? a -profiles can be physically interpreted as suitable reflecting and focusing profiles. 
From this point of view, this work can be considered as a generalization of pQ . 

Finally we note that the construction can be easily generalized to R™, n > 3. 
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